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Atomic physics was revolutionized by the development of forced evaporative cooling: it 
led directly to the observation of Bose-Einstein condensation 1 2 , quantum-degenerate Fermi 
gases 3 , and ultracold optical lattice simulations of condensed matter phenomena 4 . More re- 
cently, great progress has been made in the production of cold molecular gases 5 , whose per- 
manent electric dipole moment is expected to generate rich, novel, and controllable phases®®, 
dynamics^], and chemistry^HSl m these ultracold systems. However, while many strides 
have been made 15 in both direct cooling and cold-association techniques, evaporative cooling 
has not yet been achieved due to unfavorable elastic-to-inelastic ratios 13 and impractically 
c*5 slow thermalization rates in the available trapped species. We now report the observation 

. of microwave-forced evaporative cooling of hydroxyl (OH) molecules loaded from a Stark- 

decelerated beam into an extremely high-gradient magnetic quadrupole trap. We demon- 
strate cooling by at least an order of magnitude in temperature and three orders in phase- 
space density, limited only by the low-temperature sensitivity of our spectroscopic thermom- 
etry technique. With evaporative cooling and sufficiently large initial populations, much 
colder temperatures are possible, and even a quantum-degenerate gas of this dipolar radical 
- or anything else it can sympathetically cool - may now be in reach. 

Evaporative cooling of a thermal distributional is, in principle, very simple: by selectively 
q{ removing particles with much greater than the average total energy per particle, the temperature 

decreases. In the presence of elastic collisions, the high-energy tail is repopulated and so may 
repeatedly be selectively trimmed, allowing the removal of a great deal of energy at low cost in 
particle number. This process may be started as soon as the thermalization rate is fast enough to 
be practical and continued until its cooling power is balanced by the heating rate from inelastic 
collisions. It generally yields temperatures deep into the quantum-degenerate regime (far below 
the recoil limit of optical cooling). 



The key metric for evaporation is therefore the ratio of two timescales. The first is the rate 
of elastic collisions, which rethermalize the distribution, while the second is the rate at which 
particles are lost from the trap for reasons other than their being deliberately removed, e. g. the 
rates of inelastic scattering and background gas collisions. Both theoretica l 14 ! 17 ^ and experimental 
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wor ]$nnEM have seemed to show a generically poor value of this ratio across multiple molecular 
systems; this has led to a general belief that evaporative cooling is unfavorable in molecules^. As 
no trapped molecular system has achieved sufficiently rapid thermalization, there has been a lack 
of experiments to test this expectation. 

Hydroxyl would not, at first glance, seem to be a promising candidate for evaporative cool- 
ing. Its open-shell 2 n 3/ / 2 ground state and its propensity towards hydrogen bonding create a large 
anisotropy in the OH-OH interaction potential, which would intuitively motivate a large inelastic 
scattering rate. Chemical reactions are also possible, via the OH+OH— >H20+0 pathway; it is 
unclear whether this reaction has an activation energy barriei^lll It is thus perhaps surprising that 
the elastic collision rate actually exceeds the inelastic rate, allowing evaporative cooling. However, 
our experimental observation is unambiguous, and is further supported by quantum scattering cal- 
culations based on the long range dipole-dipole interaction between the molecules23E3 considering 
all of the fine-structure states of the rotational ground state. This analysis considers only elastic 
collision or inelastic relaxation to lower energy states, as the long-range interactions appear to fully 
dominate over short-range effects such as chemical reactions. 

In its ground state, OH has a 2 n 3 / 2 electronic character, with the lowest rotational level 
having total non-nuclear angular momentum J = \- The electronic orbital angular momentum 
couples to the rotational angular momentum to split the two opposite-parity states within J = | by a 
A-doubling of ~ 1.667 GHz; the upper parity state is labeled |/) and the lower \e). Hydroxyl is both 
paramagnetic, with a molecule- fixed moment of 2/ig (/ig is the Bohr magneton), and electrically 
polar with a dipole moment of 1.67 Debye (5.57 x 10 30 C-m). The Zeeman spectrum of OH is 
shown in Fig. Ha. Our magnetic trappy is loaded with molecules in the uppermost \ f; Mj = +§) 
state, where Mj is the laboratory projection of J. 

The results of our scattering calculations for |/; +|) molecules are shown in Fig. [Tj). The 
elastic cross section dominates the inelastic one for low energies: at a collision energy of E c = 
50 mK, the ratio of elastic over inelastic is R = 5 in a 50 mT magnetic field (B) and R = 18 in 
150 mT, while at lower energy E c = 5 mK the ratio increases to R = 23 and 137 respectively. 
Because the collisions occur in a quadrupole magnetic trap where B is inhomogeneous, the inset of 
Fig. [TJ? shows the cross sections as a function of B at E c = 50 mK. This demonstrates that inelastic 
processes are even further suppressed at B > 50 mT, in agreement with previous analysis^. 

These scattering results can be interpreted by the emergence of an effective repulsive van der 
Waals interaction between the two molecules. In zero electric field, the effective interaction in a 
scattering channel m can be evaluated in second-order perturbation theory by 

v / \ y \(m\Vdd(r)\n)\ 2 C 6 

where Vdd is the electric dipole-dipole interaction between different scattering channels, and is 
non-vanishing only between molecular states of distinct parity; and E m and E n are the asymptotic 
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energies of the relevant scattering channels. For the initial molecule-molecule channel of interest 
\f; +§) |/; +§), the other fine- structure channels are lower in energy (E n < E m ) and repel this 
energy upward. Thus, the C(, coefficient is positive, as illustrated in Fig.fTb for B = 50 mT, for the 
highest energy channel. The contribution from the next rotational state J = | of the n 3 / 2 manifold 
is too high in energy (E m — E J=5 / 2 ~ — 100 K) to give an appreciable attractive contribution at 
long-range. Hence, for low collision energies, the scattering of OH molecules is dominated by the 
long-range interaction rather than the short-range structure of the potential surface. 



Given an effective, repulsive coefficient in the incident channel, the elastic cross section 
can be approximated semi-classically at energies above the threshold regime 2 -^, as 

^ = 7t n / 5 (T(5/2)/r(3)) 2 / 5 fC 6 \ 2 ' 5 /2£ C V 1/5 (2) 



sin(^/5)r(2/5) \hj \m ied 

where m re( j = 14497 atomic units (a. u.) is the reduced mass of two OH molecules and Q ~ 
9 x 10 4 a. u. is the calculated isotropic C^. (We contrast this Q with values of 2-8 x 10 3 a. u. for 
the alkali metal atoms 27 .) Equation ([2]), plotted as a red dash-dotted line on Fig.[IJ), only slightly 
overestimates the numerical results for the elastic cross section but gives a proper trend in E c 



We observed inelastic scattering in the presence of a large electric field (Fig. [2Jd), demon- 
strating the presence of two-body collisions in our trapped sample. Motivated by the prediction 
of a favorable R, we then undertook an experiment to implement microwave-forced evaporative 
cooling. The |/; +§) — >• |e; +§) microwave transition has a small but nonzero differential Zeeman 
shift 28 , red-shifting by 26.6 kHz/mT. In the presence of a small electric field, \e; +§) molecules 
can escape the trap through the avoided crossings labeled Xj in Fig.[T^ 29 , while inelastic losses of 
\f) molecules remain unmeasurably slow. A brief microwave pulse to selectively transfer \f; +§)- 
state molecules to \e; +§), followed by a longer period with an electric field present to eject \e; +§} 
molecules from the trap, is therefore a field- (position-) selective method to remove molecules from 
the magnetic trap (see Methods). This yields what is commonly called an RF knife (Fig. |2^). A 
Zeeman depletion spectrum can be acquired by using the knife to remove molecules at a set of fixed 
frequencies and measuring the fractional depletion £; this yields a histogram of relative molecule 
number versus B. We fit this spectrum with a modified thermal distribution 

aB)dB=CoB 2 dBxcxp[--^-}x 

k B T 

fl, iffl>49.6mT 
1 exp{ — - ^ 49 f B ^ } , otherwise 

where £o is a fitting coefficient, B 2 dB oc r 2 dr is the volume element for a 3-D quadrupole trap, 
/! = 1.2 /ig is the magnetic moment of the \ f; +|) state, B is the magnetic field strength implied 
by the microwave frequency, kg is Boltzmann's constant, and T is the fitted temperature. The first 
term in the product is the simple Boltzmann expression for the molecule number as a function of 
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B, while the second is a correction for the fact that \e\ +|) molecules only disappear when they go 
through one of the X,- crossings. Specifically, while molecules transferred at fields above 49.6 mT 
(the known location of X_ 3/ / 2 ) are always energetically able to reach one of the Xj crossings and 
thus disappear, of the molecules transferred at lower fields only those with enough kinetic energy 
to fly up the trap potential to X_ 3 / 2 can escape the trap. This implies an additional Boltzmann 
factor exp[— /J.AB / hgT] in the probability of those molecules' disappearance. 

With a 5-selective technique for removing molecules, we easily implement forced evapora- 
tion by moving the edge of the RF knife from some large initial value of B towards zero, at a rate 
slow enough that the distribution remains in thermal equilibrium. Six different Zeeman spectra 
are shown in Fig. [3} demonstrating both forced evaporation and forced anri-evaporation, where for 
the latter the knife is held fixed at some low 5knife (A*#kmfe < kgT). The trap is initially loaded 
with a temperature of 51 mK; left unperturbed, it free-evaporates down to 45 mK. We have forced 
evaporative cooling by an order of magnitude down to 5.1 mK, while forced anti-evaporation can 
overcome the free evaporation and increase the temperature to 54 mK. The limit of 5.1 mK is 
attained approximately when the RF knife edge reaches X_ 3/ / 2 : while further reductions in tem- 
perature are possible, the exponential suppression of the spectroscopic signal below X_ 3/ / 2 renders 
our current system unable to measure temperatures lower than this. The plots of Fig. [3] are all 
normalized so that the area under the spectroscopic curve is proportional to the total |/)- state flu- 
orescence signal. Thus, the increase in signal height at low B in Fig. |3b-f is direct evidence of 
increasing phase- space density. 

The apparent negative signal in Fig. [3^-f can be fully fitted by assuming the presence of 
accumulated, trapped |e) -state molecules in thermal equilibrium with the visible |/) -state ones. 
Since \e; +|) molecules are totally trapped if they do not have enough kinetic energy to reach 
X_3/ 2 , they will appear, rather than disappear, during the microwave spectroscopy and contribute 
to the total depletion signal with a negative sign. Fitting the curves in this fashion gives even colder 
temperatures of 6.8 and 3.5 mK for Figs. [3^ and f, respectively. As the appearance of low-energy, 
trapped |e)- state molecules would also constitute direct evidence of evaporation, we undertook a 
direct search and indeed detected them in laser-induced fluorescence. 

We make several comments on the observed evaporation. The first is on the evaporation 
timescale: it is fast. In comparison to typical cooling rates of d{\ogT) / dt ~ 0.5 s _1 , we cool by 
an order of magnitude in only 70 ms, for a rate of 33 s _1 . This implies elastic collision rates on the 
order of 100-1000 s _1 , comparable to our trap frequenc\ 29 of ~ 1400 s _1 . As we are able to set a 
bound on the inelastic loss rate of < 2 s~ ! (as shown in Fig.[2j)), this implies a distribution-averaged 
R > 50-500, consistent with the 5-field dependence of R shown in Fig.[TJx 

We quantify the behavior of the evaporation by a set of power-law scaling relation^, three of 
which are plotted in Fig. |4j The average energy removed per molecule, rj = (trap depth) /kgT, is the 
first metric of interest: we observe a ratio r\ = 5.6, as shown in Fig. [4^. (For scale, rj =5-10 is con- 
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sidered reasonable in atomic evaporation^.) Such a relatively low value of r\ is unsurprising, given 
both the possible low (< 1000) value of the elastic-to-inelastic ratio R and the fact that molecules 
are only actually lost when they cross the spatial regions corresponding to the Xt crossings. This 
reduces the selectivity of the RF knife, as molecules transferred to the state may recollide and 
rethermalize before they find their way to a crossing. The other metrics are the behavior of tem- 
perature and relative phase-space density (PSD) as function of remaining molecule number, shown 
in Fig. [4}x The efficiency of evaporation is determined by the number of molecules remaining at 
a given temperature and PSD, that is by a = d(\ogT)/d(\ogN) and 7 = -d(logPSD)/d(logiV). 
We measure a = 1.26 and so using the fact that density scales as ^3 for our 3-D quadrupole trap 
we find 7 = 4.7. Extrapolating this 7 indicates that it would take roughly a factor of 200 reduction 
in number to increase PSD by 10 10 . 

We calculate both a and 7 assuming that our fluorescence signal is linear with total molecule 
number, justified by the optical power broadening being larger than the Zeeman broadening of the 
trap. It is very difficult to determine the sensitivity of pulsed-laser-induced fluorescence in situ, so 
we use the observed collision rate to estimate a lower bound on the density in our trap: assuming 
a unitarity-limited scattering rate /3 U with a maximum collisional angular momentum of \2h (and 
the elastic scattering rate in Fig. [TJ5 is only a factor of 3 below this value), an elastic collision 
rate of /3 u «o = 300 s _1 implies a peak density hq ps 5 x 10 10 cm~ 3 . (Integration of a Boltzmann 
distribution with this peak density over the trap suggests a total number of ~ 10 6 molecules in the 
free-evaporated trap sample and an initial peak PSD of 3 x 10~ 10 .) This density is high enough 
to permit the use of absorption imaging techniques to directly visualize the trap distribution in 
the future. Since imaging allows direct, non-spectroscopic measurement of both density and tem- 
perature, it will enable us to quantify further reductions in temperature. With the addition of a 
mechanism to prevent Majorana loss^, the favorability of R down to microkelvin temperatures 
indicates that even Bose-Einstein condensation of OH may be feasible. 

Methods Summary 

Our Stark decelerator and permanent magnet trap have been described elsewhere! 25 ! 30 ! and are illus- 
trated in Fig. [2^. Briefly, OH molecules are formed by an electric discharge through a saturated 
mixture of water vapor in 150 kPa of krypton, supersonically expanded through a pulsed valve. 
The gas packet is skimmed, focused by an electrostatic hexapole, and Stark-decelerated to 34 m/s. 
The slowed packet is then stopped at the center of the magnetic quadrupole trap by a high- voltage 
field applied between the permanent magnets, and thus loaded into the magnetic trap. 

Once trapped, the molecules are allowed to settle briefly (on the order of 5 ms) before evap- 
oration begins. The evaporation sequence consists of alternating microwave (80 /is with dBm at 
the vacuum feedthrough) and DC bias field (~650 /is at 240 V/cm) pulses: the microwave pulse 
selectively transfers hot molecules from the |/)- state to the |e)-state, while the DC bias destabilizes 
the I e) -state so that those molecules are lost from the trap. The microwave frequency is ramped 
along an exponential curve towards zero trap depth, truncated at the desired final depth. 
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After evaporation, the spectroscopy sequence is executed. It is similar to the evaporation, 
except that rather than slowly ramping a set of quasi-single-frequency pulses, each pulse is rapidly 
swept through the same narrow (50-500 kHz) frequency band with an additional sine-wave ampli- 
tude modulation (+10 dBm peak power at the feedthrough) so as to induce Adiabatic Rapid Passage 
(ARP) transference^ of all the |/) -state molecules within the frequency band to the |e)-state. A DC 
bias field again rejects the \e) -state molecules from the trap. The number of ARP pulses applied is 
generally between 5 and 30, over a 2-10 ms spectroscopy sequence. The final molecule number is 
then detected by pulsed laser-induced fluorescence using the 282 nm X 2 U^ n, v" = — > A 2 E, v' = 1 
transition. 
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Figure 1: Ground-state structure and scattering theory of OH. a, The ground- state A-doublet 
and Zeeman structure of OH. Dashed black lines are energy levels in the absence of any electric 
field. In the presence of an electric field, the level crossings become avoided (solid orange lines). 
Xi label the crossings of the \e; +§) state with the \ f; Mj = i) states; these crossings allow \e; +§) 
molecules to escape the trap in the presence of an electric field. Molecules are loaded into the 
magnetic trap in the bolded \f; +§) state, b, Cross sections as a function of collision energy: 
elastic cross section in a 50 mT magnetic field (solid), inelastic in 50 mT (dotted), and inelastic in 
150 mT (dashed). The red dash-dotted line is a semi-classical expression given by Eq. (Q. Inset: 
elastic (solid) and inelastic (dotted) cross sections as a function of magnetic field at E c = 50 mK. 
c, Adiabatic energies as a function of the inter-molecular distance at B — 50 mT. The inset zooms 
in on the repulsive van der Waals interaction for the case of two colliding |/; +§) molecules. 
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Figure 2: Trap system and inelastic collisions, a, Schematic of the Stark accelerator 30 (left) 
and magnetic trap (center rings) system. The DC-blocking capacitor (rear) decouples the high 
voltage used for trap loading from the microwave system used for spectroscopy and evaporation, 
while the laser and lens provide fluorescence detection of the trapped molecules. Inset: an artist's 
impression of evaporative cooling. As the trap depth is successively lowered by the RF knife, 
the hottest molecules escape and the remainder refhermalize to a colder temperature, b, Time- 
of-flight trace of electric field-induced collisions at 45 mK. Black squares, consistent with pure 
one-body loss, are with no applied electric field; red circles are well-fit by pure two-body inelastic 
loss induced by a large applied electric field of 3040 V/cm. (The field used in the RF knife is only 
240 V/cm, for comparison.) Error bars are one standard error. 
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Figure 3: Microwave spectra illustrating six different final temperatures. Black squares are 
data (error bars are 1 standard error), and solid lines are fits to the sensitivity-corrected Boltzmann 
form of Eq. (|3]): a the unperturbed 45 mK distribution, b anri-evaporation to 54 mK, c-f forward 
evaporation to c 18 mK, d 12 mK, e 9.8 mK, and f 5.1 mK. Dotted lines reproduce the fit from a, 
shown for comparison. Note that the x-axis scale differs between a-b and c-f; all curves have been 
shifted vertically to have a zero baseline. 
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Figure 4: Scaling relations observed in the evaporative cooling of OH. a, Final observed tem- 
perature versus effective trap depth at evaporation ramp end. b, Final temperature (solid black, left 
axis) and calculated relative phase-space density (dashed red, right axis) versus remaining relative 
molecule number. Error bars are 1 standard error; phase space density is calculated assuming a 
3-D quadrupole trap geometry, c, Molecular Boltzmann distributions implied by the fit curves of 
Fig. [3^ and c-f (right to left). Distributions are normalized so that the area under each curve is 
proportional to the total |/) -state fluorescence observed at that temperature. 
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